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Abstract 

Two equations are constructed which reflect, for fermions moving independently in a 
spherical harmonic potential, a differential virial theorem and a relation between the 
turning points of kinetic energy and particle densities. These equations are used to 
derive a differential equation for the particle density and a non-local kinetic energy 
functional. 
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1 Introduction 



There is currently considerable experimental activity in the area of harmon- 
ically confined fermion vapors [1,2]. This has been the motivation for the 
present Letter, in which equations are set up involving the particle density 
p(r) and the kinetic energy density t(r). We have written both quantities 
solely as functions of the radial distance r, as we shall be dealing with an 
arbitrary number, say (M + 1) closed shells in an isotropic external potential. 

Two equations form the basis for what constitutes a density functional theory 
of independent harmonically confined fermions in three dimensions. The first is 
a form of differential virial equation, specific to an isotropic oscillator potential 

V{r) = —muj'^r'^ (1) 
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where uP' — k/m, with k the force constant and m the fermion mass. The 
second equation relates the turning points of the particle density p(r) and 
kinetic energy density t{r), which we define from the wave function form 

t{r) ^ -^Y.^,{r)V'Ur) ■ (2) 

Two other forms of the kinetic energy density will be used below, the first 
from the gradient of the wave function 

Mr) = ^ElVV'.(r)r (3) 
and the other derived as the average of t{r) and tair), 

m = ^^^^ . (4) 



2 Differential virial equation 



As to the differential virial equation mentioned above, the early work of March 
and Young [3] , which however was restricted to one-dimensional motion along 
the X axis, say, gave by expansion of the equation of motion of the Dirac 
density matrix around its diagonal: 

This was termed a differential virial equation since, by multiplying by x and 
integrating over all x the usual integral virial theorem was recovered. We note 
from the definitions (2)- (4) above, with the consequence that 

tair) = t{r) + ^VV(r) , (6) 

that Eq. (5) can be rewritten 

dt(x) 1 , , dV(x) 

-^-^ = — p{x) — 7^ . (7) 
dx 2'^^ ' dx ^ ' 



While, for one dimension Eq. (7) is valid for a general confining potential V{x), 
its generalization of present interest refers to a three-dimensional isotropic 
harmonic oscillator with closed-shells occupancy. It is not our purpose in this 
Letter to attempt any general proof, but we shall mention later a number of 
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Fig. 1. Test of the differential virial equation: negative of dt{r)/dr (solid line) and 
{?)/2)muP'rp{r) (dots) in units fiio/af^^ as functions of the radial coordinate r/aho, 
for a 3D Fermi gas in isotropic harmonic confinement with 10, 20 and 30 filled shells 
(from bottom to top). 

convincing confirmations of our result. This is 



but with V{r) now restricted to the harmonic form (1). The factor 3/2, by 
comparison with Eq. (7), is readily understood as due to dimensionality, as 
is confirmed by multiplying Eq. (8) throughout by r to form the virial of the 
force —dV{r)/dr. A volume integration, followed by integration by parts on 
the LHS of the resulting equation yields 



where T = J d'^rt{r) is the total kinetic energy. 

It is a very simple matter to verify Eq. (8) directly for the lowest shell only 
occupied, when ipM=oi''") = ^6xp(— r^/2a^^) and PM=o{'r) = exp(— r^/a^^) 
with a/jo = ^J h/muj being the harmonic oscillator length and A a normalization 
constant. More generally, numerical tests of Eq. (8) are shown in Fig. 1 for 
three values of the number of filled shells. In these calculations we have used 
the explicit expressions for the particle and kinetic energy densities reported 
by Brack and van Zyl [4] in terms of Laguerre polynomials. 

As a last remark we note that, although the focus of the present work is on 
spherical harmonic confinement, any additional term that may enter the RHS 




(8) 




(9) 
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of Eq. (8) for more general central V{r) potentials does not contribute to the 
integral virial. Indeed Eq. (9) is the correct form of the virial theorem for such 
confinements. 



3 Turning points of peirticle and kinetic energy densities 



Having given these arguments for the form (8) of the differential virial equation 
for spherical harmonic confinement, we pass to the second central result of the 
present study, namely a relation between the turning points of kinetic energy 
and particle densities. Again, we can make a useful analogy with the one- 
dimensional harmonic oscillator. With singly occupied shells we can write, 
following Lawes and March [5]: 

the single lowest occupancy corresponding to = 1. 

We now assert that, for the three-dimensional oscillator potential energy (1) 
the appropriate generalization of Eq. (10) reads 

4^ = (M + 2);k^ - -mc^V . (11) 
p{r) 2 

Again, no general proof will be attempted. But for the lowest shell M = 
only occupied, we can again use the explicit Gaussian wave function ipM=o{f^) 
given above to calculate t(r) in Eq. (11) and p'^j^Q{r) is also readily obtained. 
Equation (11) is confirmed for M = after a short calculation. More generally. 
Fig. 2 reports numerical tests of Eq. (11) for three values of the number of 
filled shells. Again the specific expressions of Brack and van Zyl [4] for t{r) 
and p(r) have been used. 

Further confirmation, now in integral form but for general M, comes by mul- 
tiplying Eq. (11) throughout by rp'(r) and doing a volume integration. In- 
tegration by parts readily yields, invoking the equipartition theorem for the 
harmonic oscillator, 

T = -(M + 2)hw r rp'(r)4:7ir^dr . (12) 

But the integral on the RHS is related to the total number N of fermions 
occupying (M -|- 1) closed shells, the simple result 

^^l{M + 2)fuj (13) 
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Fig. 2. Test of the relation between the turning points of particle and kinetic 
energy densities: t'{r)/p'{r) in units hw as functions of the confining potential 
V{r) = for a 3D Fermi gas in isotropic harmonic confinement with 10, 20 

and 30 filled shells (from bottom to top). 

following as a consequence of Eq. (11). This is readily confirmed by calculating 
the eigenvalue sum E from the oscillator levels plus the known degeneracies. 
Again, for = 1 or equivalent ly M = 0, T = 3hw/4 follows immediately from 
the virial result T = E/2 since E = 3hLo/2 for the lowest shell. 

Before considering the consequences of Eq. (11), when combined with the dif- 
ferential virial theorem, we want to stress that, for M > 0, the Pauli Exclusion 
Principle plays a role in determining both particle and kinetic energy densities. 
Therefore, as a further confirmation of Eq. (11) we have used explicit harmonic 
oscillator wave functions (see e.g. Morse and Feshbach [6]) to construct the 
Dirac density matrix 7(r, r') for two shells occupied, i.e. for M = 1. Since the 
degeneracy of the levels for a given total quantum number n is (n+l)(n+2)/2, 
one is then dealing for singly occupied levels with 4 particles, and the Dirac 
matrix is readily constructed as 



7M=i(r, r') = {7r'/Wo)-'[l + 2{xx' + yy' + zz')/aU exp[-(r2 + r'')/2aL: 



Forming the kinetic energy density from 



(14) 



iM=i(r) = -— V^,7M=i(r,r') , (15) 

and using the explicit form of p'(r) following from differentiating with respect 
to r the diagonal density Pm=i{j') = 7M=i(r, i'')|r'=r, once again confirms the 
correctness of Eq. (11) for t'{r)/p'{r). 
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4 Density functional theory 



Given therefore Eqs. (8) and (11) as the two pillars of the present study we 
go on to construct the density functional theory of this problem of three- 
dimensional harmonic confinement. First of all, let us use Eqs. (4) and (6) to 
remove t{r) in Eq. (8) in favour of t{r) and V^p(r). Then invoking Eq. (11) the 
derivative dt{r)/dr can be eliminated, the result being a differential equation 
for the particle density p(r). This reads 



8m dr 



V^p{r) + [{M + 2)hw 



Vir)]p'ir) + ^^^p{r) = Q, (16) 



which we must again emphasize is valid for the three-dimensional potential 
energy (1). Equation (16) is evidently a third-order, linear, homogeneous equa- 
tion for the particle density p(r) of (M + 1) closed shells. Since the pioneering 
work of von Weizsacker [7], it has been clear that a major objective of density 
functional theory is to construct directly the fermion particle density from 
a given one-body potential. Equation (16) achieves that aim quite explicitly 
for closed shells in the harmonic potential (1). It is the three-dimensional 
generalization of the early result of Lawes and March [5] for one-dimensional 
harmonic confinement. 

But there is a further consequence of the two basic equations (8) and (11). 
For one can alternatively eliminate the potential, in order to directly relate 
kinetic energy and particle densities. One finds almost immediately 



dt(r) 3 , ,d 



t'jr) 
7(0 



(17) 



As emphasized earlier, t(r) and t(r) differ only by a constant times the Lapla- 
cian of the particle density. The resulting relation can be integrated, using an 
appropriate integrating factor, to give t{r) in terms of integrals of p(r) and its 
low-order derivatives. Explicitly we find 



d_ 
dr 



t'{r) 



p'[r 



)p2/3(r) 



d 



12m p^/^{r) dr 



V'p{r) . 



(18) 



One can reduce the order of the differential equation (18) by explicit integra- 
tion and by using the "boundary condition" [f (r) / p' {r)]r=o — {M + 2)fuj from 
Eq. (11) to obtain 



{M + 2)nw 
^^^^ pV3(o) P ^"^'P^"^'^ 12m 



)p\r) / 
Jo 



ds 



d 



p^/^{s) ds 



W'p{s) . 



(19) 
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A second integration leads to the final result 



t(r)^:^ + lC + ^ / ds 



3 1 3m Jo p^/^ 



Jo 



1 ^ P'is) 



2s 3p(s) 



p^/3(r) (20) 



where tw{r) = (/i^/8m)[p'(r)]^/p(r) is the von Weizsacker "surface" con- 
tribution to the kinetic energy density [7] and the constant C is given by 
C = 3(M + 2);ia;/(5p2/3(o)) - {h'^/20m){V^p{r))r=o/ P^^^{0). We have assumed 
dp{r)/dr\r=o = 0. 

The form p^/^{r) in Eq. (20) suggests, with of course the exact density p(r), 
the introduction of the Thomas-Fermi (TF) kinetic energy density defined by 

tTF{r) = cup"\r) . (21) 

The constant Ck is a multiplicative factor which we do not need to specify for 
the present purposes. Thus, as in earlier work in this Journal [8] on the one- 
dimensional case of harmonic confinement, we obtain that tTpir) and tw{r) 
enter as building blocks the non-local kinetic energy density functional (20). 



5 Summciry 



In summary, we have derived two relatively simple equations, namely Eqs. (8) 
and (11), the first being a differential form of the virial theorem and the second 
relating turning points of kinetic energy density t(r) and particle density p{r). 
These allow (i) the differential equation (16) for p(r) and (ii) the non-local 
kinetic energy density t{r) in Eq. (20) to be obtained. 
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